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On the regularity of the density of electronic wavefunctions 

S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and T. 

0stergaard S0rensen 


Abstract. We prove that the electronic density of atomic and molecular 
eigenfunctions is smooth away from the nuclei. The result is proved with¬ 
out decay assumptions on the eigenfunctions. 


1. Introduction and main results 

We consider an JV-electron molecule with L fixed nuclei whose non-relativistic 
Hamiltonian is given by 


N 

(1.1) H = %(R,Z)=^ 

.7=1 


L 


Zi 


3 \ x i ~ Rl \, 


+ E 



E 

1 <l<k<L 


ZiZ k 

I Ri-RkY 


where R = (R±, R 2 , ..., Rl) € R 3L , Ri ^ R k for k ^ l, denote the positions 
of the L nuclei whose positive charges are given by Z = (Zi, Z 2 ,. ■ ■, Zl). The 
positions of the N electrons are denoted by {x\, X 2 , ■ ■ ■, Xn) € R 3Ar , where Xj 
denotes the position of the j-th electron in K 3 and A = A 

dimensional Laplacian. 

For shortness, we will sometimes write (0 as 


H = -A + V. 


The operator H is selfadjoint on L 2 (M. 3N ) with operator domain TF 2,2 (R 3Ar ). We 
consider the eigenvalue problem 

(1.2) H = Ei/j, 

with if) £ L 2 (K. 3Ar ) and E the real eigenvalue. We will assume, without loss, that 
the eigenfunction if} is real valued. For i£l 3 define 

■ ■ . ,Xj-i,X,Xj+i, .. .,x N ), 
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and 

d£j — dx i • • • dxj_idxj + i ■ ■ ■ dxpf. 

Then the density associated with ip is defined as: 

N 

(1.3) 

j=i ' 

Theorem 1.1. Let ip satisfy & Then the associated density p satisfies 

(1.4) P G C°°{R 3 \{R u ...,R l }). 


N f 

p(x )= sy ^ d± -j- 


Remark 1.2. We remark that analogous natural smoothness results hold for 
the 2-electron density and the 1-electron density matrix, see Q. 


Remark 1.3. The electron density plays the central role in various approxima¬ 
tion schemes for the full 31V-dimensional Schrddinger equation for IV-electron atoms 
and molecules (Thomas-Fermi, Hartree-Fock, and Density Functional Theory). 

Strangely enough, the density p has only rarely been the subject of mathe¬ 
matical analysis, see [li], Q, and references therein. Our present results are not 
surprising—surprising is the fact that it has not been shown some decades ago. 
Since the eigenfunction ip is real analytic away from the singularities of the poten¬ 
tial one would even expect p also to be real analytic away from the nuclei. 


Remark 1.4. Theorem 1.1 is an improvement of Jy Theorem 1] since we have 
no decay assumptions on ip. This improvement will come about by using Sobolev 
spaces instead of Holder spaces. 


For simplicity we will here only explicitly consider the case of atoms, i.e. we 
take L = 1 and can without loss of generality define R\ = 0. To prove Theorem 0 
it clearly suffices to prove 

(1.5) pG<7°°(R 3 \{0}) 

with p redefined by p(x) = J ip 2 {x, X 2 , ■ ■ ■, Xn) dx 2 • • ■ dxw, which will be done 
below. 

Let us recall the definition of Sobolev spaces (with integer exponent). 

Definition 1.5. For fceN and an open set D C K" the Sobolev space W k ’ 2 ( Ll ) 
is defined as 

W k ’ 2 { n ) = {/ G L 2 (D)| d a f G L 2 {n) for all a G N", |a| < k}. 


2. Differentiation of ip ‘along’ singularities of V 

The key to the smoothness of the density is the understanding that one is 
allowed to differentiate ‘a long ’ singularities of the potential. The precise statement 
of this is given in Lemma |2.2| . The statement and proof of this lemma is analogous 
to a similar result in |j]. Proposition 2], the difference being that we here measure 
regularity of functions by demanding that they be in W 2,2 , whereas the analogous 
space in [|lj. Proposition 2] was the Lipschitz continuous functions C 0,1 . 

Definition 2.1. Let D C 1" be an open set. We define R°°(D) as follows 
= {/ G C7°°(R n )| supp/ C Q and d a f G L°°(R”) for all a G N n }. 
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Lemma 2.2. Let P,Q be a partition of {1,..., N} satisfying 
P/0, P fl Q = 0 PUQ = {1,...,N} 


Define, for P , Q as above and e > 0 


( 2 . 1 ) 


Up(e) = {(aq,.. .,x N ) G K 3W 


\xj | > e for j G P, 

Xj - x k \ > e /or j G P, k G q|. 


Define also 


x P 


1 

vW 


jep 


and let T be any orthogonal transformation of R 3Ar smc/i that T{x\,... ,Xn) = 
{xp,x') with x' G R 3iV ~ 3 . Let ip satisfy ([h^) and let <p G B°°{Up{e)). Then 

d* p ({<jnp) o T*) G W 2 ’ 2 (R 3Ar ) for all a G N 3 . 


We may assume without loss of generality that P = {l,...,iVi}, with Ni < N. 
Then the orthogonal transformation T can be written as: 

f 1 ... 1 n • • • 0 

( 2 . 2 ) T= l VNl ,/Nl 

\ 1 1 • • • tpj 1 . tN 


with the first row being understood as 3 x 3 matrices—first N\ repetitions of -^=/ 3 
and then N—N\ repetitions of the 3x3 0-matrix. The remaining part of the matrix, 

(2.3) f = ( ti ■ ■ ■ tN ) with tj G M 3 jv_ 3i3 (R), 


is such that the complete matrix T is orthogonal. 

Let us write ip = ip o T*. Notice that, since A is invariant under orthogonal 
transformations, ip satisfies 


(2.4) 


-Aijj + (V o T*)ip = Eip. 


We can now prove Lemma 2.2. 


Proof. We will prove that for all </> G B°°(Up(e( 1 — 2 fc ))) we have df p {(pip) G 
W 2 , 2 (R 3Ar ) for all 7 G N 3 with |y| < k, with <j = <p o T*. The proof is by induction 
with respect to k. For k = 0 there is nothing to prove since D{H) = W 2 , 2 (K 3Ar ). 
Let k > 0 and suppose that df. p {(pip) G W 2,2 (K 3Ar ) for all <p G B oc {Up{e{ 1 — 2~ k ))) 
and all 7 G N 3 with I 7 I < k. Take now 7 G N 3 with I 7 I = k + 1. Let <p G 
B co {Up{e{ 1 — 2 _ P" +1 )))). From the eigenfunction equation (2.4) for ip we get 

-A (4>ip) = (-2 vipvip - 4>A4> + Eijyip) 

N 


(2.5) 


-E 

3= 1 


-Z 


of 


(H>) - E 

1 <j<k<N 


1 


\Xj - x k \ 


o T* I {(pip). 


Now we differentiate equation (|T^) with respect to xp (in the distributional sense). 
By the induction hypothesis it is clear that 

dl p {-2\'4>Vip-ipA4> + E4>ip) e L 2 (R 3N ). 


Notice that d xp commutes with \xj\ 1 o T* for j G Q and with | Xj — Xk\ 1 o T* 
for j,k G P or j,k G Q. Furthermore, on supp</> the remaining parts of the 
potential, |xj | _1 o T* for j G P and Xj — Xk | -1 o T* with j G P, k G Q, are 
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bounded functions with bounded derivatives of arbitrary order, so we get, using 
the induction hypothesis that 

dl p {(VoT*)W) = (VoT*)dl p (W)+g, 

where g £ TT 2,2 (R 3Ar ). Furthermore, by assumption we have d^ p {ftft) £ IV 1,2 (R 3Ar ). 
Therefore we obtain (V o T*)(82 p (ftft)) £ L 2 (R 3N ) since the operators and 

|Xj — £fc| _1 are bounded from W 1,2 (R 3Ar ) to L 2 (R 3N ) (see for instance [3[ p. 169]). 
Thus, in the sense of distributions, ( |2.5| ) implies 

A 82 P (ftft) G L 2 (R 3N ). 

Via standard elliptic regularity results, this implies that 82 P (ftft ’) £ W 2,2 (R 3Ar ) and 
finishes the proof. □ 


3. The differentiation of p 

In this section we perform the differentiation of the density, as stated in (El- 
Let ft £ Cq°(R 3 \ {0}). We will prove that 

(3.1) 82(ft 2 p) £ L 1 (R 3 ) for all 7 e N 3 . 

By a Sobolev embedding theorem, we therefore get (ft 2 p £ C'°°(R 3 ), which implies 
( 0 ), since ft £ Cq°(R 3 \ {0}) was arbitrary. 

We use the partition of unity introduced in [jl]. Let R > 0 be such that 
{x £ R 3 | jo) < Rjdsuppft = 0. Letxi 5 X 2 be a partition of unity in R + : X 1 +X 2 = 1, 
with Xi(a;) = 1 on [0, f?/(4/V)], suppxi C [0, R/(2N)) and Xj G C , °°(R + ) for 
j = 1,2. We combine the yys to make a partition of unity in R 3W . Obviously, 

i= n (xi(i* 3 --* fc i)+x2(i*i-* fc i)). 

1 <j<k<N 

Multiplying out the above product, we get sums of products of Xi’s and X 2 ’s. We 
introduce the following index sets to control these sums: Define first 

M = {(j,k)£{l,...,N} 2 \j<k}, 

and let 

/CM, J = M\I. 

Now define, for each pair /, J as above, 

fti( x )=( II Xi(kj-a*l))( II X 2 (K-ar fc |)). 

(j,k)ei U,k)eJ 

Then we get 

1 = II ( Xl + X2 )(\ x j ~ x k \) = Y ^/( x )> 

1 <j<k<N /CM 

where the sum is over all subsets I C M. 

Therefore we have, with gj = ft 2 ft 2 , 

p{x 1 ) = ^2 9i(xi,x 2 , ... ,xn) d,X 2 • • • dx N = ^ pi(x i). 

ICM ^ /CM 

To verify (0 we have to prove that for every / C M, 82(ft 2 pi) £ L 1 (R 3 ) 
for all 7 £ N 3 . As in |lj, for each I C M we can choose a P (depending on 
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I), with 1 £ P, and such that supp </>(xi)</>/(xi,... , xjv) C Up(R/4N) (using that 
supp^fl{ |x| < R} = 0). This will be essential in the considerations below in order 
to apply Lemma [2.2[ 

Let u £ Cq° (IF) . Choose / £ Co°(M 3Ar_3 ), with / = 1 for |x| < 1, and let 
/k(x) = f(x/K), K > 0. We calculate 


J R 3 
(3.2) 


(d% u)(xi)(4?pi)(xi) dxi = / (d2 1 u)(x 1 )((/) 2 g I ){x 1 ,...,x N )dx 1 ---dx N 

Jr 3N 


= ,7|/2 / ( d2 p u)(4> 2 gi){x P ,x')dxpdx 

Jm? N 


accor ding to the coordinate transform given in Lemma |2.2| as specified in (2.2) 
and (2.3). Here u = u o T* and sim ilarly for the other functions. By Lebesgue 
integration theory it follows that (3.2) equals 


(3.3) 


N~ hl/2 lim 

K —KX) 


R3JV 


( h K 4> 2 g I )(x P ,x') dxpdx ', 


with h K {x P ,x') = ( d2 p u)(x P ,x')f K (x 

Noting that xp = y/Nl{x\ — t\x'), t\ being the transposed of t\. and that u and 
Jk are Cq° functions, hx £ C^°(R 3JV ). Hence (3.3) equals 

(3.4) 


(—l)l' y l J /V 1 ' 7 ^ 2 lim f ufxv dxpdx ', 

K^CO J R 3N 


with v(xp,x') = d’]. p {(p 2 gi){xp,x'). Here the derivative is taken in the distribu¬ 
tional sense. 

Suppose we have shown that v £ L 1 (R 3Ar ), then by Lebesgue integration the¬ 
ory (3.4) equals 


(3.5) 


(— 1 ) 1^1 N ' 7 ^ 2 f uv dxpdx ', 

J r 3N 


and by transforming again and applying Fubini’s theorem (3.5) equals 

(3.6) (_i)ItIj\t“ItI / 2 f f LW dx 1 , 

Jr 3 

wi th w ( x i) = f R 3 N -3 (v o T)(xi, ..., xjv) dx 2 • • • dxx, w £ L 1 (R 3 ). Combining (3.2) 
(3.6), (|3.l| ) follows. 

So to finish the proof of dll) it remains to verify that v £ L 1 (R 3JV ), which will 
follow via Lemma 2.2. In order to apply the lemma we need the following. 


Proposition 3.1. Let u,v £ L 2 (R n ) and d Xl u, d Xl v £ L 2 (I 
sense). Then 

d Xl {uv) = (d Xl u)v + ud xi v, in the distributional sense. 


(distributional 


Proposition 3.1 is easy to prove by standard density arguments. 

Remember that v = d) p (</> 2 4> 2 ip 2 ) , and that P was chosen such that supp (jxf>i C 
Up(e) for e = R/4N. Let now v\ = Then v = d) p (vf). Since £ 

B°°(M 3N ) it follows via Lemma 2.2 and Leibniz’ rule that d xp (v\) £ L 2 ( R 3Ar ) for 
all a £ N 3 . From this and Proposition [O] it is easy to prove by induction that 

E (^a P ®l)(^a P ®r) £ i 1 (R 3Ar ). 

a+/3=7 

This finishes the proof of (|l^) and hence of Theorem 1.1 for the atomic case. 
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